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Abstract 



We compute the first-order correction to the correlation functions of the 
stationary state of a stochastically forced harmonic chain out of equilibrium 
when a small on-site anharmonic potential is added. This is achieved by 
deriving a suitable formula for the covariance matrix of the invariant state. 
We find that the first-order correction of the heat current does not depend 
on the size of the system. Second, the temperature profile is linear when the 
harmonic part of the on-site potential is zero. The sign of the gradient of the 
profile, however, is opposite to the sign of the temperature difference of the 
«_! ■ two heat baths. 



1 Introduction 

The goal of this paper is to begin a perturbative analysis of invariant probability 
measures arising in the context of non-equilibrium statistical mechanics. As a model 
at hand, we will consider a Hamiltonian chain of N oscillators interacting through 

1 Supported by the Japan Society for the Promotion of Science and by Les services federaux des 
affaires scientifiques, techniques et culturelles du gouvernement Beige . 
2 Supported by the Swiss National Science Foundation. 
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nearest-neighbour interactions, coupled at its boundaries to stochastic heat baths 
of different temperatures, and that we will perturb by a small anharmonic (quartic) 
on-site interaction. The covariance of the stationary state in the purely harmonic 
case has been computed in ^21 QUI- For other cases, i.e. anharmonic cases, almost 
nothing is known about the physical content of the stationary state, except results 
about the positivity of entropy production and validity of linear response theory [U] . 

It is a natural idea to attempt to understand its physical properties by performing 
a perturbative analysis. Such an approach, based on the phonon picture, has been 
exploited by physicists to tackle the Fourier law, see pQ for a classical exposition. In 
particular, the Peierls theory seems successful in computing the thermal conductivity 
and its thermal and dimensional dependence. The Peierls approach assumes from the 
beginning the existence of an infinite non-equilibrium state where local temperature 
equilibrium is expected to hold. It is also based on several implicit assumptions, 
such as the validity of a Boltzmann equation for phonons. In this paper, we adopt 
a different approach and begin a rigorous perturbative analysis of a finite (although 
taking N large will have some simplifying features) anharmonic chain. Our starting 
point is a formula, which seems to be new, for the correlation functions of the 
stationary measure. This formula allows us to derive (matrix) equations for the 
first-order correction. The relationship between our approach by stationary non- 
equilibrium states (SNS) and the Peierls approach is, at this stage, far from clear. 
A first interesting step would be to achieve some understanding of the equivalence 
of the definition of the thermal conductivity by the Green-Kubo formula and its 
definition in the SNS approach as, roughly speaking, the ratio of the heat current 
and the temperature gradient. 

The main obstacle to developing a perturbative expansion of SNS's is that, in 
contrast to the equilibrium case, no explicit formula for the invariant density is 
known. Moreover, the fact that the relevant models are degenerate in a stochastic 
sense makes it laborious to obtain a systematic perturbative expansion starting from 
the equations of motion. We circumvent this difficulty by deriving a formula for the 
two-point correlation functions of invariant states, which holds under the assumption 
of /^-convergence of the finite-time correlation functions to those of the (unique) 
invariant measure. We emphasize that the validity of the formula is not restricted 
to the concrete problem of the anharmonic chain considered here. It may prove 
useful whenever the invariant measure is not explicitly known, in particular in the 
context of transport phenomena modeled by hypoelliptic stochastic processes. We 
also remark that the form of the formula for the covariance is very similar to, and 
provides a lower bound on, the expectation of the Malliavin matrix. 

Our main result concerning the heat current is that its first-order correction re- 
mains uniformly bounded as the number of oscillators goes to infinity. In particular, 
perturbative analysis does not, at first order, reveal any sign that Fourier law holds 
in such anharmonic models as numerical studies suggest, see e.g. jH]. Furthermore, 
we find that the first-order correction to the temperature profile is exponentially de- 
caying in the bulk of the chain, with a decay rate that depends on the strength of the 
harmonic part of the on-site potential. When this strength vanishes, the correction 
to the temperature profile is linear. However, the sign is "wrong" , in the sense that 
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the linear profile has the lowest temperature near the hottest bath and the highest 
temperature near the coldest bath. This is analogous to the result of [12], where the 
temperature profile is also oriented in the "wrong" direction. The main difference is 
of course that in , the temperature profile is exponentially decaying. In order to 
understand what is responsible for this awkward behaviour, it would be interesting 
to examine the perturbation theory of harmonic chains that are defined differently 
near their ends, e.g., with respect to the harmonic interaction or the coupling with 
heat baths. Another feature of our solution is that the temperature profile is shifted 
downwards, in the sense that the temperature at the middle point of the chain is 
lower than the arithmetic mean of the temperatures of the heat baths. 

The remainder of this paper is organized as follows. In Section 2, we specify 
the basic set-up for the type of anharmonic chains we will consider. Section 3 is 
devoted to the derivation of our basic formula for the covariance. In Section 4, we 
derive the matrix equations for the first-order corrections to the harmonic case. One 
assumption of this section is that the invariant measure is regular in the anharmonic 
parameter. We postpone the proof of this fact to a future publication. The last two 
sections are devoted to the resolution of these equations. This is done by generalizing 
the methods of HOI H21 - 



2 A model for heat conduction 

In order to explain the behaviour of the thermal conductivity in crystalline solids, 
one often models the solid by a chain (or lattice in higher dimension) whose ends are 
coupled to heat baths maintained at different temperatures. The coupling can be 
taken stochastic and more precisely of Langevin type. In one dimension, the set-up 
is as follows. At each site of a lattice {1, . . . , N} is attached a particle of momentum 
Pi and position g^. The dynamics is Hamiltonian in the bulk and stochastic through 
the Langevin coupling to heat baths at the boundaries. The Hamiltonian is of the 
form, 

N N 

H (hl) = E(jP? + V ^)) + E Ufa - + U( qi ) + U(q N ). (2.1) 

i=l i=2 

Specific choices for the potentials U and V will be specified below. The equations 
of motions are given by, 

dq t = Pi dt, % = 1, . . . , N, (2.2) 
dH 

d Pi = -—{p,q)dt, i = 2,...,N-l, (2.3) 
oqi ~ ~ 



and, 



dpi = —Q^-(p,g)dt - -ypidt + yj2jkT t dw t , (2.4) 
dH 



dpN = — q — (p,q)dt — ^p N dt -\- J2^kT N dw r . (2.5) 
oqN 
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T\ and Tn stand for the temperature of the left and right reservoirs, respectively, 
whereas wi and w r are two independent standard Wiener processes. 

It is an easy fact to check that when 7\ = T/v = T = the measure on 
the configuration space R 27V whose density with respect to the Lebesgue measure is 
given by 

p(p,q) = Z- 1 e- pH ^D (2.6) 

is invariant (stationary) for the stochastic dynamics defined above. Explicitly, one 
can check that for L the generator of the dynamics and any function / in its domain, 

J Lfp(p,q)dpdq = 0. (2.7) 

In the case of two different temperatures, existence, uniqueness and exponential 
convergence to an unique invariant state has been established under fairly general 
conditions on the potentials U and V El HI] • In the case of harmonic coupling, 
the covariance of the stationary state has been exactly computed in fT2\ HUj. 

An essential ingredient of the proof of the uniqueness is the fact that the system 
satisfies the so-called Hormander condition. This condition implies that the noise 
spreads in a sufficiently good way through the system, so that the transition proba- 
bilities have smooth densities. This property is encapsulated in the non-degeneracy 
of the Malliavin matrix associated to the stochastic system under study. As the noise 
represents the injection of energy into the system, it is natural to enquire about the 
relationship between the Malliavin matrix and the correlation functions of the sta- 
tionary state. This might provide a way to tackle the description of the stationary 
state when its density is not explicitly known. Indeed, from a physical point of 
view, the central question, once uniqueness has been established, is to compute the 
energy spectrum and correlation functions of the stationary state and ultimately, to 
establish the validity of the Fourier law. As mentioned above, the case of a harmonic 
chain has been completely and explicitly solved. The main feature of the solution is 
a flat temperature profile and an associated infinite thermal conductivity. 

The basic idea in order to perform a perturbation theory of the non-equilibrium 
stationary state is to write the two-point correlation function of the stationary mea- 
sure under a "Malliavin" form, similar to the form derived by Nakazawa in the 
Gaussian harmonic case, |lUj . 

3 The Malliavin matrix and the covariance ma- 
trix of the stationary measure 

We consider now a general system of stochastic equations. Denote by x t E R d the 
solution of the stochastic differential equation, 

n 

dx t = X (x t ) x k{xt) dw k {t) (3.1) 

k=l 

with initial condition xq = x, where the toj-'s are n independent one-dimensional 
Brownian motions and Xi, I — 0, ... ,n, are C°° vector fields over R d satisfying for 
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any mult i- index a, 

||a a XKa;)||<C(l + ||x||^) (3.2) 

for some K a > 0. We note that solutions to such equations are in general not ensured 
to exist globally. In the sequel, we restrict ourselves to the following situations. 

Assumption 3.1. For all x G R d , equation 113. 1\) has a unique strong solution xt, 
t > 0. This solution has finite moments of all order: for all p > 1, T < oo, and 
x G R d , there exists a constant C = C(x,p, T) < oo such that for < t < T, 

E x (\\x t \\r) <C. (3.3) 

When in need of emphasizing the dependence of the solution to (|3.1|) on the initial 
condition x and the realization of the d-dimensional Brownian motion w in the 
interval [0, t), we shall write it as x t (x,w([0,t})). We denote by V t the associated 
semigroup, 

V t f(x)=E x (f(x t )) = J f(x t {x,w([0,t})))^(w([0,t)), (3.4) 

where P is the <i-dimensional Wiener measure, by A the generator of the semigroup, 
and by L the associated second order differential operator, 

d d 

L = J2*odi+Yl ^ijdidj, (3.5) 

i=l *J=1 

where, with ® denoting the tensor product, 

n 

a= 1 -Y / X k ®X k . (3.6) 

fc=i 

From Assumption 3.1 on the process solution xt and the bounds ()3.2j) for the vector 
fields X[, it follows that for each t and w[0, t], the map x h- > x t (x, w[0, t}) is C°° on R d 
with derivatives of all orders satisfying the stochastic differential equation obtained 
from (|3.1|) by formal differentiation. Furthermore, for all multi-index a, p > 1, and 
t > 0, 

E{\\&*x t (x,-)\\ p ) <oo. (3.7) 

In the sequel, we will denote U t (x,w[0,t\) = Dx t (x, w[0, t}), where DX denotes the 
Jacobian matrix of a vector field X on R^. The matrix Ut is the linearized flow and 
it solves the equation, with initial condition Uq = 1, 

n 

dU t = DX (x t )Ut dt + J2 DX k (x t )U t dw k (t) . (3.8) 

k=l 

Below, E x U t denotes fU t (x,w[0,t])dP(w[0,t]). 

Let us now assume the existence of an invariant probability measure \i for the 
process solution x t of (|3.1|) and consider the covariance matrix at time t, 

C t (x) =E x (x t ®x t ) -E x x t ®E x x t . (3.9) 
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The following result is the starting point of the perturbative analysis performed in 
subsequent sections. It provides an expression for fi(Ct) in terms of the linearized 
flow Ut, where is a shorthand notation for j Rd /(x) dfj,(x). 

Proposition 3.2 Suppose that the bounds A3.2\) and Assumption 3.1 are satisfied. 
Suppose in addition that the invariant measure \i for the process solution x t of A3.1]) 
is such that the functions x \— > E x x\ , x i— > LE x x!,, and x i— > ciij(x)E x U 3 s l , belong to 
L 2 (R d , do) for all i,j, I, and s <t. Then, 

fi(C t )= / dsJ2KVU s X k (.)(g>EU s X k (.)). (3.10) 

Proof. We will show below that the map s i— > /i(Ei s ® E x s ) is different iable, with 

d n 

— //(E.x s ® Ex.) = - ^ ^(E. £/ a X fc (.) (8) E £/ s X fc (.)) . (3.11) 
ds k=i 

Identity (|3.1U|) thus follows from the invariance of the measure //, since 

ti(C t ) = fi{E.{x t ® X*)) - |i(Ei t <g> E xt) (3.12) 
= ^(x <g) x) - /i(E x t <g) E x f ) (3.13) 
/•* d 

= -/ rfs— /i(Ex s ® Ex s ). (3.14) 
Jo ds 

To obtain (|3.11|) . we first note that (j3.3j) implies that any function / G C 2 (R d ) 
with first derivatives of at most polynomial growth is in the domain of the gener- 
ator A with Af = Lf. Similarly one easily checks that for such /, (|3.7j) implies 
•A(V t f) = L(Vtf). Therefore, Kolmogorov equation yields -^(E x x s ®E x x s ) = 
LE X x s Cg) E x x s + E x x s £g> LE X x s , which, by Holder inequality and our assumptions, 
belongs to L 1 (R a! , dfi). Thus, 

^(E.i s <g> Ex s ) = u.(LE x s <g> Ex s + Ex s <g> LEx s ). (3.15) 

Let us next define for /, g G C 2 (R d ), 

T{f,g) = L(fg)-fLg-gLf, (3.16) 

which reads 

T(f,g) = 2 y £ «,,<),/ <),<!. (3.17) 

Since it follows from ()3.7|) that c^E^x^ = E x lip, our assumptions imply as above that 
r(E.x l s , Ei J s ) G L 1 (R d , rf/i) for all z, j. It follows in particular that L(E.x s ®E x 8 ) G 
L 1 (R d , d/i) . Because of the invariance of /i (which implies fi{Lf) = 0), we are thus 
free to subtract from the /^-expectation on the right hand side of 1)3.15)1 a term 
L(Ex s ® Ex s ), so that 

^((E.x s <8> E,x s )ij) = -/i(r(E.x l s , E x J j). (3.18) 
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Formula ()3.11j) finally follows from the computation, recalling (jH.fij) . 

n 

V{K.y.Kri)(-r) = J2 (^xU s X k (x) ® B x U s X k {x)) . . . (3.19) 

This concludes the proof of Proposition 3.2. 
Proposition 3.2 immediately implies the 

Corollary 3.3. Suppose that the hypothesis of Proposition 3.2 are satisfied for all 
t > 0. Suppose in addition that 

lim C t = nix ®x)- ixix) <g> nix) = $ , (3.20) 

t — >oo 

in L 1 (R d , d/j,). Then, 

$ = / dsY,v(VUsX k (.)®EU s X k (.)). (3.21) 
Jo k=i 

The expression f!3.21|) for the covariance matrix of a stationary state is the basic 
formula that we shall use to develop a perturbation expansion in the next section. 
Since both sides of (|3.21j) involve an averaging with respect to /i, it is not clear at 
first sight how informations on /i can be extracted from ()3.21|) . We observe, however, 
that in the case of a linear drift Xq and constant vector fields X k , k = 1, . . . , n, all 
expectations may be dropped and (|3.21|) becomes 

$ linear = f" ds U s (T X k ® X k ) Uj. (3.22) 

J ° k=l 

One thus recovers the standard formula for the covariance of the stationary state of 
a linear stochastic equation with constant diffusion coefficients. As we shall see in 
the next section, it is possible to iterate this simple observation in order to begin a 
perturbation expansion. 

Another feature of formula (|3.1U|) is to provide a link between the covariance 
matrix C t and the so-called Malliavin matrix. The Malliavin matrix associated to 
equation ()3.1|) at time t reads, in the normalization of 0, 

f t n 

M t = dsJ2 U t V s X k {x s ) ® U t V s X k {x s ) , (3.23) 
Jo k=i 

where V s is the inverse matrix of U s . An easy computation reveals that fi(EJM t ) can 
be expressed in a form closely related to (J3.1Uj) . namely, 

n(EMt) = dsY, v(E.(U s X k (.) <g> U s X k (.))). (3.24) 
Jo kTi 

Indeed, we first observe that for s > fixed, Y* = U t V s satisfies Y s s = 1 and 

n 

dY s l = DX (xt)Yt dt+J2 DX k {xt)Y2 dw k {t) (3.25) 

k=l 
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for t > s. Comparing with ()3.8|) yields that Y" s * = Yg(x s (x,w[0, s]),w[s,t]) has the 
same P-distributions as U t - S (x s (x, w[0, s]), w[s, t]), where w{r) = w{r) — w(s) for 
r > s. Furthermore, for x fixed the map w i— > Yg(x,w[s,t]) is w[0, s]-independent. 
Therefore, since (x,w) i— > F s '(x, w)Xfc(x) <E> w)Xfc(x) is measurable, one may 
use the Markov property of Xf to write, 

E x (F s *(x s )X fc (a; s ) ® r s *(x s )X fc (x s )) = E se (E I ^ e .(£/ t _ a ( 3 /)X fc ( 3/ ) ® U t - B {y)X 1i {y))). 

(3.26) 

Identity (|3.24p then follows by using the invariance of the measure /x and changing 
variables in the integral over s in ()3.23|) . As a consequence, Proposition 3.2 provides 
a lower bound on the expectation of the Malliavin matrix. 3 

Corollary 3.4. One has 

(i(C t ) < fi(EM t ). (3.27) 
Proof. The inequality simply follows from (|3.1U|) . (|3.24j) . and the matrix 



U s X k {x) - E x U s X k (x)) ® (C7 a X fc (x) - E x U s X k (x))] (3.28 



being positive definite. 



4 Perturbative analysis of the non-equilibrium an- 
harmonic chain 

We shall analyze the effect of adding an anharmonic perturbation to a modification 
of the model treated by Rieder, Lebowitz and Lieb [12]. We consider the case of a 
harmonic chain with fixed ends to which one adds an anharmonic on-site potential, 
i.e. in f)2.1jl . we set 

U{x)= 1 -u 2 x 2 and V = \u 2 kx 2 + ^Ax 4 . (4.1) 



The model considered in [12] has k = but the computation of the covariance of the 
stationary state is very similar and the result is given below. We write the equations 
of motions ()2.2|) - (J2.5j) under the matrix form, 

with J\f(q) and dw the vectors in H N given by J\fi(q) = qf and dwj = 5uy/2jkTi dwi+ 
5N i ^2 r ykT N dw r , and 

H-L-a) 

3 The order relation is defined in the following way. For two matrices Xi,X2, we say that 
Xi > Xi whenever X\ — Xi is a positive definite matrix. 
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where g K and a are N x N matrices given by (g K )y = ^ 2 ((2 + «)<% — 5 

ij+l Oij—l) 

and 3Lij = jSij(S\j + $Nj)- Above, 1 denotes the unit matrix and the zero matrix 
or vector, as is clear from the context. We note that the stochastic terms in (|4.2j) 
are given by constant vector fields, namely, in the notation of Section 3, 

Xk= (d k ) where (dfc)j = <W 2 ^ T *> ( 4 - 4 ) 

for k — 1,N. In particular, the coefficients involved in the generator L are 
constant. They are given by 

E X k ®X k =( ° ° ) , (4.5) 
fc=i,jv V / 

where A^- = 2 r yk5 i j(T 1 5ij + T n 5nj). Furthermore, the linearized flow U x of (|4.2|) is 
given by 

dU x = bU x dt - 3\C x (t)U x dt , (4.6) 

where 

CHt) = ( J (t) I ) , (4.7) 

with v X j(t) = 5ijqf(t) and qi(t) the gj-component of the solution of (|4.2|) at time t. 
Finally, we note that the matrix b in (14. 2J) has the property that all its eigenvalues 
have strictly negative real part. A proof of this fact can be found in JU] modulo 
obvious modifications. 

In order to study perturbatively the SNS of our chain, we would like to use the 
identity (|3.21|) . However, some of the hypothesis of Corollary 3.3 related to the 
invariant measure are not known to hold for equation ()4.2|) when A > 0. (The case 
A = has been covered in |12j.) Although from a mathematical point of view, this 
is not a mere technical problem, but since the main goal of this paper is to illustrate 
the use of formula (J3.21)) for perturbative analysis on a specific example, we will 
assume that these hypothesis hold, see Assumption 4.1 below and the remark that 
follows. On the other hand, Assumption 3.1, i.e., the existence of strong solutions 
and their moments, follows from standard techniques and we briefly discuss it now. 
We first note that for A > 0, the function H(q,p) = 2N + H(q,p), with H the 
Hamiltonian given by (|2.1|) and ()4.1|) . satisfies 

i%,p)>C(l + ||g|| 2 + |H| 2 ), (4.8) 

for some C > and all (q,p) G H. 2N . Thus, H is a C 2 (R 2N ) confining function. 
Furthermore, one computes 

(LH)(q,p) = - 7 (p 2 + p%) + 2 7 A;(T 1 + T N ), (4.9) 

which implies that LH is uniformly bounded by above. A classical result, see e.g. [7], 
Thm 4.1, then ensures for all initial conditions (q,p) G R 2JV the existence of a unique 
global strong solution to ()4.2|) . Regarding the bounds (j3.3J) . they are an immediate 
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consequence of the following a priori bound. For any 9 < (2k max{Ti, T/v}) , one 
has 



E (g,p) 



< e 2 7 fc0(T 1 +r iv )i e eH(q, E ) _ (4-10) 



Bound (j4.10j) can be obtained in a similar way as in the proof of Lemma 3.5 in jTTj . 
However, the existence of a unique invariant measure for (|4.2|) is still an open prob- 
lem. We thus introduce the following 

Assumption 4.1. The finite time truncated two-point correlation function of the 
process defined by A4.2\) converges to the covariance matrix of a unique stationary 
measure /i A in L l (H 2N , dfi x )-norm. Furthermore, the decay properties of /i A are such 
that E (£i p)[( M )], LE (M) [(q t ,p t )}, and E {g ^[U x ] belong to L 2 (R 2N , d/i x ). 

Remark. The uniqueness of the invariant measure is proved in (HI E] for a large 
class of anharmonic chains. The invariant measure has a smooth density with ex- 
ponential decay and is shown to be mixing 4 . An important restriction is that the 
potential U must not grow asymptotically slower than V, and thus equation (|4.2jl 
does not fall into the class covered in 0[H]. However, as is argued in [TT], the fact 
that the on-site potential grows faster than the nearest-neighbour interaction should 
not affect the ergodic properties of the measure but only the rate of convergence. 
Although we could consider a similar anharmonic chain with an additional quar- 
tic term in the nearest-neighbour interaction, the equations that one then needs 
to solve, see below, are computationally more involved. Furthermore, restricting 
to (|4.2j) will allow us to compare our results to the usual A0 4 expansion when the 
temperatures of the two baths are equal. 

Provided Assumption 4.1 holds, let $ A denote the covariance matrix of the unique 
stationary state of equation (|4.2jl and express it according to (|3.21j) as 

$ A = / dt ]T v x (EU x X k <g> BU x X k ) . (4.11) 

Jo k=l,N 

We first briefly review the harmonic case A = 0. As mentioned at the end of the 
previous section, one obtains from (j4.11|) 

$° = / dte ht T>e h \ (4.12) 

where 



o 



D= £ X k ®X k = ( ° ° j , (4.13) 

k=l,N \ J 

with Ay = 2jk5ij(Tx5ij + T^5Nj)- Since the eigenvalues of b have strictly negative 
real part, the integral in (|4.12jl is convergent and it follows from integrating by parts 
in b$° that $° must satisfy the equation 



b$° + $°b T = -D. (4.14) 



4 In the result is actually stronger. The convergence to the unique invariant measure is 
shown to be exponential. 
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The unique solution of this equation has been explicitly derived in [12]. It is given 
by 

$° = 



X z 



where, denoting T = I)J f K 1 V = and G K = w~ 2 g 



and 



<S>° y = kT(l + v Y°), 

$° = — r/Z°, 

7 



/ 01 

02 
03 



V o 



>N-2 <PN-1 



\ 



>N-l 



-01 / 



(4.15) 

(4.16) 
(4.17) 
(4.18) 



(4.19) 



o 



(4.20) 



/ 

-01 



V 



W-2 0JV-1 \ 
07V-2 



•W-l 



<?1 

/ 



Above, v = ^2 and the quantities 0j, 1 < j < N — 1, satisfy the equation 

JV-l 

i^v+n )ij4>j = $U , 

3=1 



(4.21) 



(4.22) 



where denotes the fc-square matrix given by (G^l K )ij = (2 + is + n)5ij —5i ; j + i — 

fiij-l- The solution of (|4.22|) is given by 

siixh^-iV 
J sinhiVa 

with a defined by cosh a = 1 + [v + k)/2. Hence, one has for large N and fixed j 
the asymptotic formula 0j = e _QJ . In the context of SNS, one usually defines the 
temperature to be the average kinetic energy, i.e. in our case, 



(4.24) 
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It is easy to see that the above solution yields an exponentially flat profile in the 
bulk of the chain. 

We now turn to the first-order perturbation of the anharmonic chain. We first 
introduce our second assumption on the process solution of (|4.2j) . 

Assumption A2. The measure fi x is absolutely continuous with respect to the 
Lebesgue measure and as a function of A its density p x (x) is C°° in a neighbourhood 
of 0. For all x, all derivatives are bounded in a neighbourhood of 0. 

Remark. The proof of this fact should follow from an analysis similar to the ones 
developed in [3] or to prove the smoothness of the probability transitions in a 
parameter of the related stochastic differential equations. 

To derive an expression for $ x = ^$ a |a=o ; we compute from (j4.11j) 

$ x = JU A |a=o (4.25) 



i=l,N 
roo 



+ dt J2 E— f/ t A | A=0 X,(.)®E.^ X 4 (.))+tr., (4.26) 



i=l,N 



and observe that the first term vanishes because /i 1 = ^/i A |A=o integrates constants 
to zero. In order to compute the last terms, we first evaluate W t = ■^■U x \\=o . 
Deriving with respect to A on both sides of equation ()4.6|) . we get 

dW t = bW t dt-3 C°(t) U? dt, (4.27) 

from which it follows that, since Wo = 0, 

W t = -3 C ds e h{t - s) C°{s) e hs . (4.28) 
Jo 

Inserting (J4.28j) in (|4.26p . we obtain, using in addition the invariance of 

roo rt 

$! = -3 / dt ds V e h{t - s) Ne hs Xi®e bt Xi + tr., (4.29) 
Jo Jo i=hN 

= -3 r dt [ ds e b( *" s) N e hs D e bT * + tr. , (4.30) 
Jo Jo 

where D is given by (|4.13|) and 



Exchanging the integrations over t and s and changing variables leads to 

$ x = -3y rfteb ' N (y o ^e bs De bTs )e bTl + tr. , (4.32) 
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which, with (|4.12j) . finally yields, 



$! = -3 / dte ht (N$° + $°N T )e b * . (4.33) 
Jo 

The method used to derive the above equation will also provide the equations for 
the next orders of the perturbative expansion. However, obtaining them concretely 
requires some more work and we reserve that part and the general Feynman rules 
for a further publication. We note that integrating by parts in (J4.33)) yields the 
equation for $ x 

M) 1 + $ 1 b T = 3(N$° + $°N T ). (4.34) 

In Section 6, we will derive an explicit expression for $ x and thus for the first order 
correction to the heat current and temperature profile. It turns out to be easier to 
do so by solving equation (J4.34|) rather than by using (J4.33j) . In the next section, 
we first make a few preliminary remarks about equations of the form (J4.34D . 



5 Solving the equation for the first order 

The symmetry properties of the inhomogeneous term in equation 1)4.34)1 will play a 
special role. We will need to consider symmetry properties both with respect to the 
diagonal and to the cross- diagonal. 

Notation. For a i^-square matrix M, we denote by M c the transpose of M with 
respect to the cross- diagonal, namely, (M )^- = M# +i~j 1 K+i-i ■ 

Definition. We call a square matrix M c-symmetric or c- antisymmetric if M c = M 
or, respectively M c = — M. Denoting 



1 

1 



(5.1) 



we call a 2N -square matrix M CT-symmetric or CT- antisymmetric if M c = JMJ 
or, respectively, M c = —JMJ. 

We first list a few properties of equations of the form ()4.34j) . 

Lemma 5.1 Let b as above and H a 2N-square matrix. One has: 

(a) . The unique solution of the equation 

b$ + $b T = H (5.2) 

is given by 

$ = -f(Iie b 'He bTt . (5.3) 
J o 

(b) . If H is CT-symmetric or CT-antisymmetric, then $ is CT-symmetric or, re- 

spectively, CT- antisymmetric. 
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(c). If H is of the form 



then the solution of (15.2)) is of the form 



H= ?! , (5.4) 



•= - X Z Y V < 5 - 5 > 



Proof. Point (a) follows from the matrix b having all its eigenvalues with strictly 
negative real part. Indeed, this property implies that the operator $ i— > b$ + $b T 
is invertible, and integrating by part in b<3> reveals that (|5.3|) is the unique solution 
of (|5.2|) . Point (c) is obvious, whereas (b) follows from the identity Jb c J = b T and 
uniqueness of the solution of ()5.2|) . 

Lemma 5.1 implies in particular that $ 1 is the unique solution of ()4.34J) and is 
of the form 

*' - ( -k t\ ) ■ ^ 

In particular, it follows from (|5.fi|) and $ 1 being symmetric that is antisymmetric. 
In order to find an expression for the solution of equation ()4.34j) . we decompose the 
inhomogeneous term on the RHS of ()4.34|) into powers of rj and solve the equation 
separately for each case. One has 

3PT 2 

3(N$° + $°N T ) = :L ^(H + r^Hi + r? 2 H 2 ), (5.7) 



where, cf. (OSJ - lOEt and IjOTt . 



H o-(v o °G- G 'o V °)' ^ 
1 " lv ViG- 1 + V X° 7 z.[V , Z°] J ' (5 - 9j 
H2= (v X 7 z/jvrZ ])' (5 ' 10) 



with 



V = diag(G; 1 ), V 1 =diag(X°). (5.11) 
In the sequel, we will denote (Vo)y = S^g^ where g>j = (G,: 1 )^ read 

sinh ia sinh(iV + 1 — i)a /K <riN 
smha smh(iv + l)a 

with a defined by cosh a = 1 + k/2. Writing 

3A- 2 T 2 

9 1 = -^-{91 + ^1 + ^1), (5-13) 



CJ 4 
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one thus obtains that $j }, I — 0, 1, 2, is the unique solution of 

b®} + <^b T = B_i . (5.14) 
In order to scale out the constants in b, we denote for I = 0,1,2, 

= ( i x z , f, ) ■ 

together with 

R = 7 " 1 a, G K =^ 2 g K , (5.16) 

namely, R^ = 6ij(8ij + S^j) and (G re )y = (2 + k)^- — 5ij + \ — 8ij-\. The zero order 
term in (J5.13|) is just the first-order perturbation of the anharmonic chain at the 
equilibrium 7\ = T^. Inserting (J5.15|) into (|5.14|) for I = yields the equivalent 
system of equations for X , Y and Z 

Y = X G K + Z R + G; 1 V , (5.17) 
[G K ,Z ] = -^{R,Y }, (5.18) 

with the requirement that Xo,Yo are symmetric and Zq is antisymmetric. One 
easily checks that its unique solution is given by 

X = -G^VoG^ 1 , Y = 0, Z = 0, (5.19) 

thus recovering, as expected, the first-order correction of the A0 4 model. Proceeding 
similarly for $J and $^ one finds that Xi, Y 1; Zi solve 

Yx = XxG K + ZxR + (X°V + G- 1 ^), (5.20) 
[G K , Zx] = -i{R, Yx} + [Z°, V ], (5.21) 

whereas X 2 , Y 2 , Z 2 solve 

Y 2 = X 2 G K + Z 2 R + X°V 1 , (5.22) 
[G K , Z 2 ] = -i{R, Y 2 } + [Z°, Vx]. (5.23) 

Furthermore, using the c-symmetry properties of the solution X° and Z° of the 
harmonic case, cf. (J4.19|) and (J4.21j) . one easily checks that Hi is CT-antisymmetric, 
whereas H 2 is CT-symmetric. This implies that Xx, Yx are c-antisymmetric and Z x 
is c-symmetric, whereas X 2 ,Y 2 are c-symmetric and Z 2 is c-antisymmetric. This 
simply reflects the fact that changing the sign of rj corresponds to interchanging the 
reservoirs at the ends of the chain. 

In the next section, we will derive explicit expressions for the solutions of the 
above equations. To this end, we will need the following identities. Let X be a 
solution of 

[G K ,X\=U, (5.24) 
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with U a given matrix. It thus follows from [G K , X]y = Uij that 

X^j+i — Xj_x,j = Z^y + (Xj + i ) j — Xj j3 -_x), (5.25) 

where matrix elements with an index equals to zero or N + 1 are set to zero. Let 
us first consider X antisymmetric. In particular, X is entirely determined by its 
elements Xjj with i < j and satisfies Xj+i^ — X^j-i = — (Xjj+i — Xj_x,j)- For 
i < j, applying ()5.25j) recursively j — i times thus leads to 

1 

X u+ i - Xi_i )3 - = - ^U i+ ij-i . (5.26) 
z z=o 

This gives all matrix elements Xy, 1 < j < N. Applying ()5.26|) recursively z — 1 
times finally leads to 

= 9 X! X! Mi+i-k,j-i-k-i ? (5-27) 

for z, j such that i < j. Proceeding similarly, one obtains for a c-antisymmetric 
matrix X satisfying ([5.24)1 . 

j i-l N-i-j 

= n E X! Mi+l-k,j+l+k+l , (5.28) 

for i+j<N. If X is both antisymmetric and c-antisymmetric, one iterates identity 
f!5.26j) N + 1 — i — j times to obtain 

1 J— *— 1 N-i-j 

Xjj = — - X! X] Mi+l+k+l,j+l-k , (5.29) 
4 fc=0 z=o 

for i < j and i + j < N. Finally, proceeding similarly but without assuming any 
symmetry properties, one derives an expression for X depending both on U and the 
first line of X, 

i i—l i—k 

Xjj = 2J X lji+ j_2fc+l — ^^W i+ i_ fc _ u _ fc+i , (5.30) 
fe=l fc=l Z=l 

for 1 < z < j and z + j < N + 1 . Formula (|5.3U|) will be used later for X symmetric 
and c-symmetric. It reflects the fact that in such cases, the solution of (|5.24j) is 
determined up to a polynomial P(G), that is up to iV independent variables which 
can be supplemented as the first line of X. 

6 The first-order correction 

In this section, we derive an expression for the first-order correction to the heat 
current and temperature profile. We find that the part corresponding to the heat 
current is uniformly bounded in N. In particular, a first-order perturbation does not 
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reveal any sign that Fourier law might hold in such anharmonic models, as numerical 
studies indicate, see e.g. jH]. Indeed, if Fourier law holds whenever A is finite, one 
might expect the derivatives of the heat current to develop a singularity at A = 
when iV — > oo. 

Regarding the temperature profile, the part of the solution proportional to rj is 
exponentially decaying in the bulk of the chain whenever k > 0. The decay rate is 
slower than in the purely harmonic case. For k — 0, the profile proportional to 77 
is linear in the bulk of the chain and we compute its slope explicitly. However as 
explained in the introduction, the sign is "wrong" , in the sense that the linear profile 
has the lowest temperature close to the hottest bath and the highest temperature 
close to the coldest bath. The same type of phenomenon is present for k > 0, see 
Figure 1. Moreover, we observe that the part proportional to i] 2 gives a significant 
contribution, which results in a shift of the temperature at the middle point of 
the chain. The temperature at this point is no more the arithmetic mean of the 
baths temperatures. Although surprising, this is a phenomenon which seems to be 
observed in numerical studies of certain anharmonic chains, see [8 . 



6.1 First-order correction to the heat current 

In our model, the heat current in the SNS is given by ($3)^+1. The first-order 
correction will thus be given in terms of, cf. ()5.13|) and (j5.15|) . 

3£* 2 T 2 

*i = j-CZo + ijZi + ifZa). (6.1) 

By ()5.19|) . Z does not contribute and one easily checks that for 1 < i < N — 1, 

(Z 2 ) M+1 = 0. (6.2) 

That is, Z2 does not contribute to the current either. Indeed, recall that Z2 is 
antisymmetric and satisfies equation ()5.23j) . Since {R, Y 2 } is a bordered matrix 
and [Z°, Vi] is zero on the diagonal, one obtains by using formula (J5.27|) that 

(Y 2 )n = (Z 2 )i2 = (Z 2 )23 = • • • = (Z 2 )./v-ijv • (6-3) 

v 

On the other hand, the c-antisymmetry of Z 2 implies that (Z 2 )i 2 = — (Z2)jv-ijv, 
which leads to ()6.2j) . We note for later use that this also implies 

(Y 2 ) u = 0. (6.4) 

It thus remains to consider the contribution of Zx- Since Zi is antisymmetric, one 
obtains from (|5.21j) that 

Z 1 = Z + Z, (6.5) 

where Z and Z are given by formula ()5.27|) with U replaced by — ^{R, Yi} and, 
respectively, [Z°, V ]. We first observe that {R, Yi} is a bordered symmetric matrix, 
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so that formula (|5.27jl yields 

/ o 

-<pi 



(fN-2 <PN-1 \ 



V>1 



(6.6) 



-<P2 



-<P1 



J 



where the quantities ipi, . . . ,<fN-i are related to the first line of Y 1; namely, for 
j = l,...,iV-l, 

i/^ = -(Y 1 )y. (6.7) 

Furthermore, [Z°, Vq] having zero diagonal implies that -Z^i+i = 0. One therefore 
obtains 

(Zi)i,i+i = Z M+1 = tpi . (6.8) 

In order to compute the vector (p e R^ -1 , one considers the first line of equation 
(|5.2(Jj) for Yi into which one substitutes identity (|6.7|) . We first need to compute 
Xi. Equation ()5.20|) and the symmetry properties of Xi, Yi and Zi imply that Xi 
satisfies 



[G K , XJ = {R, Z x } + ([X°, V ] + [G- 1 , Vx]) 

= {R, Z} + {R, Z} + ([X°, V ] + [G-\ VJ). 

Since Xi is c-antisymmetric, it follows from (j6.10p that 

Xi = X + X, 



(6.9) 
(6.10) 



(6.11) 



where X and X are given by formula ()5.28j) with U replaced by {R, Z} and, respec- 
tively, {~R,Z} + ([X°, V ] + [G~\ Vi]). Using that {R, Z} is a bordered antisym- 
metric matrix, one obtains from (J5.28j) and (J6.6)) that 



¥2 



X 



V?2 



V -VN-1 



(fiN-2 fN-l \ 

— (PN-l 



~<P2 
V?2 ~fl ) 



(6.12) 



Equation (J5.2U)) now reads 



with 



Yi = XG K + ZR + W, 
W = XG K + ZR + (X°V + G^Vi), 



(6.13) 
(6.14) 
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and since (XG K + ZR) y = (G^XiOj = (G^V); for j = 1, . . . , N- 1, where G^ 
denotes the /c-square version of G K , it follows from (|fi.7j) that 

G^V=-w, (6.15) 

where w e R^ -1 is given by Wj = Wy, j = 1, . . . , AT — 1. Therefore, one finally 
obtains, recalling that r\ — Ti ~ Tn 



2T ' 



(*JW = V , (6-16) 

with 9? given by = — [Gy+ K 1 ^] _1 w. As represent the first-order correction 

to the current, it is consistent to see that they are all equal to each other. 

Before turning to the first-order correction of the temperature profile, we study 
the behaviour of <pi with N. We first note that X solves the equation [G K , X] = 
{R, Z}, as is easily checked from ()6.6|) and (|6.12jl . This implies that X solves, 
cf. (16~TU|) and (I6~TT| . 

[G K ,X] = {R,Z} + ([X°,V ] + [G~\ Vx]), (6.17) 

which in turn implies, by using in addition the symmetry properties of the matrices 
involved in (|6.14|) . that W is c-antisymmetric and satisfies the equation 

[G K , W] = G K ZR + RZG K + (G K X°V° - V°X°G re ). (6.18) 

Hence, Win = and it follows from formula ()5.28j) that 

w = w « + W ( 2) , (6.19) 



where, for 1 < j ' < N — 1, 



w ■ 



(2) 



1 N-j 

- "£(G K ZR + RZG K ) lJ+j , (6.20) 
1 i=i 

- (G K X°V° - V°X°G K ) M+J - . (6.21) 
1 i=i 



We first consider w^ 1 ). We note that G K ZR + RZG K is a bordered c-symmetric 
matrix and that Z is c-symmetric since both Z\ and Z are c-symmetric. One thus 
obtains from ()6.20|) 

w (i) = QVf-i)^ (6 _ 22) 

where, for 1 < j < N — 1, 

2, = Z ljj+1 . (6.23) 

In order to compute Z, we note that Z solves the equation [G K , Z] = — ^{R, Yi}, 
as is easily checked from (|6.6j) and (|6.7|) . Therefore, Z solves, cf. (|5.21|) and (|6.5j) . 

[G K ,Z] = [Z°,V ], (6.24) 
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and since Z is antisymmetric, as both Zi and Z are, it follows from (j4.21|) . (Vo)y = 
Sijgi, and formula (|5.27J) . that for 2 < j < N, 

13-1 

Z ^ = 9 - ^)0j-2z , (6.25) 

with the convention cf)_k — — <j>ki < k < N — 1. Thus, is given by ()6.22|) with 
Z G R^ -1 given by 

1 i 

Z i = o J2(9j+i-i ~ 9i)<f>j+i-2i • (6.26) 
z i=i 

We next consider w^ 2 \ We first note that 

G K X°V - V X°G K = (G U+K X°V - V X°G V+K ) + u(V X° - X°V ), (6.27) 
and compute, using ()4.19|) . ()4.22|h and (V )y = S^gi, that for i < j, 

(Gh.«X V - VoX^^)^- = Sugrfj-i + 8 Nj g i( f> N ^. (6.28) 

Therefore, 

(G K X°V - V X°G K ) ii = 8iigj4>j-i + 8jN9i<t>N-i + v{g% - gj)<pi+j-i , (6.29) 

with the convention 4>N+k = —<pN-k, < k < N. One thus finally obtains for 
w^ 2 ) G R^" 1 , using in addition that gn-j = 9j+i, 

v N ~ j 

w j = 9j+i<Pj + ^J2(9i~ 9j+i)^j-i+2i ■ (6.30) 
A i=i 

Using (JETS), (l6~TTHl . jS23, (|S2SJ, l|Oty . and the fact that the 0/s decay exponen- 
tially, it is easy to see that is uniformly bounded in N. 

6.2 First-order correction to the temperature profile 

We now analyze the first-order correction to the temperature profile. It is given by 
($J)« where, cf. (EZBD and j£B), 

= — j-(Y„ + 17Y1 + r/ 2 Y 2 ). (6.31) 

By (|5.19p . Y does not contribute to In order to compute the diagonal of Y 1; we 
use the fact that Yi is c-antisymmetric and satisfies the equation, as a consequence 

of J52DJ, 

[G K , Yx] = G.ZxR + RZxG. + (G K X°V - V X°G K ). (6.32) 

Using (|5.28|) . (|6.29j) . and the fact that g 2 i = gN-2i+i, one thus obtains for 1 < i < 
[N/2], where [x] denotes the largest integer smaller or equal to x, 

N-i i-1 



(Yi)« = (Gf -^ZOa-! + (g 2i 2i _! + - £> 2 z £(<fc-* " 0+*fi)), ( 6 -33) 
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9 

- ' fc=0 



where Zi G R w 1 is given by (Zi)j = (Zi)ij+i. Since the 0j decay exponentially fast 
with rate a, see ()4.23j) . it follows that all terms but the first give an exponentially 
flat contribution to (Yi)u. We thus write, and will adopt a similar notation in the 

(Y 1 ) u = (Gf -^Zi)*-! + 0(e- a J). (6.34) 

In order to compute the dominant term in the above expression, we first use that 
Zi = (f + Z where Z is given by f)6.26|) . and G^+^V = — w where w = Q^~^Z + 
w (2) with w (2) given by ijOJt . to obtain %_ = (G^+«^) _1 (i/^ - w (2) ) and thus 

(Yi)« = ((G^V'Gf - X) (^ - wW))^ + 0(e-*). (6.35) 

It follows from the expression (l6~2Tl) for w( 2 ) and properties of G^" 1 ), G^ 1} , and 
their inverse, that the second term gives an exponentially flat contribution to the 
temperature profile. To compute the remaining term y = ^(G^~ 1 ^) -1 G^ V ~ 1 ^ Z , we 
first note that it satisfies 

G^y = ^Gf- 1 ) Z. (6.36) 

We next compute G^~^Z. In the expression (|6.26j) for Z, changing the summation 
index to k with 2k = j + 1 — 21 if j is odd and 2k = j — 21 if j is even, one obtains, 
using in addition the symmetry properties of g i: that for j > 2 

^ , E fc ii(^i±i +fc - g^- k )<p2k if j is odd, 3 
ELi(ffi+fc - ^i+i-fc)^- 1 if J is even - 



2 ' 2 



For j = 1, Z\ — 0. Computing the differences of g's arising in the above expression 
leads to 

~ _sinh(iV-j)a * S mh(2k - j)a 

Z > ~ sinh(iV + l)a ti sinha 02fe " J ' (6 ' 38) 
where J = if j is odd and J = 1 if j is even. Hence, Z can be rewritten as 

Z, =n s ^fJ^ + (e-">), (6.39) 
3 n sinh(JV + l)a v ; 

where the constants po and pi are given by 

'g? Btoh(»- g) a L 

smh a 

A straightforward computation finally leads to, recalling that cosh a = 1 + k/2, 

(Gf -^), = (2 + K )( Pl - Po) s S ^~j^ + C x 5 Xj + <D(e-<v), (6.41) 

where C\ is a constant that depends on N and a only. It thus remains to compute 
the vector y given by equation ()6.36|) . To this end, we note that a vector of the 
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form ()6.41j) is almost an eigenvector of G[,+ re 1 ' ) . More precisely, one has for v with 
Vj = (-lp"+ 1 sinh(A^- j)a, 

(Gl+,; 1) f) i = (4 + v + 2k) Vj + Sjj sinh Na. (6.42) 

Therefore, writing 

J (4 + z/ + 2 K ) sinh(iV+l)« + J ' 1 j 

and inserting in ()6.36|) yield for r the equation (G^^r)j = C 2 5ij + C(e _aj ) with 
C*2 a constant depending on N and a, cf. (j6.41|) and ()6.42|) . whose solution reads, 



by using (JQ2]) . 

rj = C 2 ^ + 0(e- Qi ). (6.44) 

Hence, r is an exponentially decaying correction to y as given by f)6.43j) . Finally, 
since (Yi)^ = yn-\ for 1 < 2 < [N/2], we obtain from ()6.43|) . 

, v n K 2 + «)(Pi -po)sinh(iV+ 1 -2i)a _ 2ai 

(Yl)tt = (4 + , + 2.) sinh(iV + l)a + ° {t } ' (6 ' 45) 

Since Yj is c-antisymmetry, ()6.45|) also gives the elements (Yi)jj for [iV/2] + 1 < 
% < N. In particular, since cosh a = 1 + k/2, it follows that the contribution of Y x 
to the temperature profile is exponentially flat in the bulk of the chain whenever 
k > 0. When k = 0, on the other hand, a = and Yi gives a linear profile. In the 
limit N — > oo, it is straightforward to compute that for k = 0, pi and po are given 
by 

Po = 2 — anc ^ Pi = 2 — ' (6.46) 

2 sinh a 2 sinh a 

with a defined by cosha = 1 + u/2. One thus has pi — Po = 1/(4 + v) and the 
temperature profile for k = is given by 

CYx)« = 77T~~T9 (lw~7T - l)+0(e" 2Qi )- (6-47) 



(4 + z/) 2 VjV + l 

The temperature profile is linear, but oriented in the "wrong" direction. Indeed, if 
for instance T\ > Tjy, then one obtains from (|6.3ip . which involves a multiplication 
by rj = (T\ — T/v)/(Ti + T/v), that the slope is positive. 

We next consider the contribution of Y 2 to the temperature profile. Since Y 2 is 
c-symmetric, it will introduce, if nonzero, a global shift in the temperature profile. 
As we shall see, this is indeed the case. To compute the diagonal (Y 2 )a, we proceed 
as for Yi. We first recall that (Y 2 )n = 0, cf. ()6.4|) . and note that Y 2 also satisfies, 

[G K , Y 2 ] = G K Z 2 R + RZ 2 G K + (G^X ^ - V 1 X°G K ). (6.48) 

Denoting by r/> the first line of Y 2 , i.e., 

V»i = (Y 2 )ii, (6.49) 
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1 20 40 60 80 100 



Figure 1: Contribution of Yi to the temperature profile (is = 1, N = 100). 

one uses (15.30(1 to obtain from (|6.48(1 the following expression, for % > 2 and 2i < 

N + l, 

i—l i—l k 

(Y 2 )« = ^2fc+i - Uk-i+i,k+i , (6.50) 

k=l k=l 1=1 

where ipx = (Y 2 )n = has been used, and 

U = G K Z 2 R + RZ 2 G K + (G K X V! - V!X G K ). (6.51) 

Since Y 2 is c-symmetric, (|6.50|) determines all diagonal elements (Y 2 )jj, 2 < % < 
N—l. The first term on the RHS of ()6.48j) is a bordered matrix and a straightforward 
computation yields 

k 

XXGkZ 2 R + RZ 2 G K )fc_/ + i i fe + ; = (Q K C)2k i (6.52) 
i=i 

where £ denotes the first line of Z 2 , i.e., 

Ci = (Za)w- (6-53) 

The second term on the RHS of (|6.51(l is identical to the corresponding term ap- 
pearing in (|6.18() . with V replaced by the diagonal matrix (Vi)ij = 5y02i-i- For 
1 < i < j < N, it is thus given by, cf. (IFT2H . 

(G K X°Vi - ViX?G K )ij = v(^-i - fcj-i^i+j-i 

+ Sa4>2j-i4>j-i + 5 jN 4>2i-i4>N-i , (6.54) 
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with the convention (pN+k = —<pN-k, < k < N. Inserting (|fi.52j) and (|fi.54j) into 
(FKHjl leads to 

(Y 2 ) u = £ A fe (6.55) 

k=l 

where, for k > 1 and 2k < N — 1, 

Afc = ^ 2 fc+l - (G K C)2fc - (02fc-104fc-l + ^2fc ^2(<f>2(k-l)+l ~ 4>2(k+l)-l)) ■ (6.56) 

1=1 

One checks that |Afc| decays exponentially. First, recalling (|4.23|) and our convention 
(pN+k = —<t>N-ki < k < N, this is clearly true of the last two terms in (|6.56|) . Next, 
an expression for the first line of Y 2 can be obtained from equation ()5.23j) by using 
that Z 2 is c-antisymmetric. Formula ()5.28|) and (Z 2 )k,k+i = 0, cf. ()6.2|) . imply that 
for 1 < k < [(iV-l)/2], 

^ ^ k N-k-1 

-1p2k+l = -Z {<h(l+n)+l ~ 02(J-n)+l) , (6.57) 

V 1 n=l l=k 

with the convention 4>N+k = —(pN-k, < k < N. In particular, ^fe+i decays 
exponentially. We finally compute (, the first line of Z 2 . One has Ci = Gv = by 
antisymmetry and c- antisymmetry of Z 2 , and applying formula ()5.29j) to equation 
(EH) yields for 2 < j < N - 1 

I i-i N-j 

= 7 E ^i-2n Z (#2('+n)-l - 02(j+i-n)-l) , (6.58) 
4 n=l i=l 

with the conventions (fi-k = —(pk and 0Ar + fc = — <pN-k, < k < N. Therefore, one 
has for 2 < i < [(JV + l)/2], 

(Y 2 ) ll = /i + 0(e- ai ), (6.59) 
where the constant h is given by 

h = h l + vh 2 , (6.60) 

with 



h = Z ( 2 C2fc+l - (2 + K)( 2 k - 02fc-104fc-l), (6.61) 

fe=l 

F— 2~ 1 k 

Z (-^fc+i - ^2fc 53(02(*-J)+i ~~ 02(fc+Q-i)) • (6-62) 



fe=l >zy i=l 

A straightforward, but lengthy, computation yields the following asymptotic formu- 
las for large N, 

cosha(cosha — 1 — k/2) 

hl = o a • P ' ^ 6 ' 63 ) 

2e a sinn a smh 3a 

1 / 1 cosh a \ 

^2 = 5— ( — ; — + ; )• 6.64 

4sinh a v cosh a e a smh 3a 7 
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Figure 2: Contribution of Y 2 to the temperature profile {y = 1, N = 100). 



Recalling that cosh a = 1 + iy + k)/2, one obtains 



h 



(u + k)(2 + v + k)(4 + z/ + k) ' 



(6.65) 
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